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The exact normalization of a multicomponent generalization of the ground state 
wavefunction of the Calogero-Sutherland model is conjectured. This result is ob- 
tained from a conjectured generalization of Selberg's A^-dimensional extension of the 
Euler beta integral, written as a trigonometric integral. A new proof of the Selberg 
integral is given, and the method is used to provide a proof of the multicomponent 
generalization in a special two-component case. 



1. INTRODUCTION 

The quantum many body system (Calogero-Sutherland model) is the subject of 
much present day interest due to its connection with quantum chaos [1] and fractional 
statistics [2,3]. The development of these applications has been greatly assisted by the 
recent discovery [4,3] of mathematical methods which provide the exact evaluation of 
ground state correlations, both static and dynamic. These exact calculations rely on A^- 
dimensional integration formulas, which are generalizations of Selberg's [5] A^-dimensional 
extension of the beta integral: 
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\l=l''^ J l<j<k<N 



r(Ai + 1 + jA)r(A2 + 1 + jA)r(i + jX) 

/=o r(Ai + A2 + 2 + (Ar + j-l)A)r(l + A) ^ ■ ' 

Generalizations of the Calogero-Sutherland Hamiltonian to include internal degrees of 
freedom of the particles have recently been formulated [6]. These models are of interest 
in condensed matter physics because of their relationship with quantum lattice models, 
notably the exchange t — J Hamiltonian [7] . However, from the viewpoint of exact 
calculations, the theoretical development of these models is not as advanced as that of 
the original model. In particular, the analogue of the Selberg integral (1.1) for the exact 
multicomponent ground state wavefunction [6] 
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is not known in the existing htcraturc. 

This deficiency has motivated us to pursue the task of formulating the appropriate 
analogue of the Selberg integral. We first transform the Selberg integral into an equivalent 
form due to Morris [8]: 
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where CT denotes the constant term in the Laurent expansion. The constant term is well 
defined when the function is a Laurent polynomial: this requires a + &, c G Z>q, a.b E Z. 
However the Fourier integral in (2.3) defines an analytic function for Rc(a+6) > 1, Rc(c) > 
(at least) and is evaluated by the product of gamma functions whenever it is defined. 
By using a combination of analytical and numerical methods we are able to conjecture 
the value of (1.3) when the product of differences is replaced by (1.2), and the one-body 
factors extended over all particle coordinates. We are able to provide analytic checks 
on the validity of our conjecture. In fact, in the case p = 1,A = 1 these checks suffice 
to prove the conjecture. Our results include the conjectured exact normalization of the 
wavef unction (1.2). They can also be cast in a form which directly generalizes the Selberg 
integral (1.1). Furthermore, they also include the conjectured exact normalization of the 
ground state wavefunction 
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for the system in a harmonic well. 



2. EQUIVALENCE BETWEEN SELBERG-TYPE INTEGRALS AND 
SOME FOURIER INTEGRALS 
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2.1 The inter-relationship 

To perform a numerical investigation of the value of the Selberg integral (1.1) or any gen- 
eralizations, it is convenient to first transform the Selberg-type integral into an equivalent 
Fourier integral. For this purpose we use the following lemma. 

Lemma 1 

Let f{ti, . . . , tN] {p}) be a Laurent polynomial in ii, . . . , In, with {p} as parameters. For 
Re(e) large enough so that the r.h.s. exists 




This result follows immediately from term-by-term integration of the Laurent polyno- 
mial. Note that for e an integer the Fourier integral is equal to 

N 

CT{t„...,t^}l[nfi-ti,...,-tN;{p}) (2.2) 

1=1 

From this lemma we see that the Selberg integral with Ai arbitrary and A2, A non- 
negative integers is equivalent to the Fourier integral in Morris's integral (1.3) with a — b 
arbitrary and a + b,c non-negative integers. 

2.2 Numerical evaluation 

For e an integer the Fourier integral in (2.1) can be computed by exact numerical inte- 
gration. Thus, whenever 



Pi Pn 

2iTi{xini+...XNnN) 

"ni,...,niv^ 

m=-pi nN=-PN 
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we have 

/ dxi... dxNg{xi,...,XN) = j-j- ...j-j- J2 dini/Mi^.-.^UN/MN) (2.3) 



provided Mi > pi {I — 1,...,N). This result follows by term- by-term integration and 
summation of the Fourier series for g{x) . In both cases only the rij — term {j — 1, . . . , N) 
of g{x) remains. 

We will use the formula (2.3) below to provide numerical data on the evaluation of 
some generalizations of (1.3). 

2.3 Analytic properties of the Fourier integrals 

Denote the Fourier integral on the l.h.s. of (2.1) by /(e; {p}). Suppose furthermore that 
/ is real when each 9j is real so that 

/(-e^-'^S . . . , -e^-'^^; {p}) = f{-e-^^''\. . . , -e'^-^^^; {p}) (2.4) 

Some immediate properties of /(e; {p}) are 
(i) /(e; {p}) is an entire function of e. 
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(ii) J(e; {p}) = /(— e; {p}) (this property requires (2.4)). 

(iii) I{e;{p}) = ((sin7re)/7rj J{e;{p}), where J(e;{p}), which is given by the r.h.s. of 
(2.1) for Re(e) large enough, is a rational function of e. 

Consider further property (iii). Suppose in fact that J(e; {p}) is the reciprocal of a 
polynomial: 

J{^; M) = mux: ^^^^^^ (2-5) 

Then qj{{p}) is the maximum number of times the integer nj{{p}) occurs as a power in 
a single term of the Laurent expansion of /. We will see below that the Selberg-type 
integrals related to the wavefunction (1.2) have the property (2.5). 

In the cases that (2.5) holds, a conjectured evaluation of (2.1) of the form 

((sin7re)/7r) a({p})Ajv(e; {p}) (2.6) 

can be proved to be correct up to the multiplicative function of the parameters Q;({p}) 
by an application of Liouville's theorem. Thus we need to show that 1/A]^{t] {p}) is a 
polynomial in e and to specify the positions and orders of its zeros. Then we need to show 
that the positions of the poles of JAr(e; {p}) coincide with these zeros, and their order is 
no greater than that of the corresponding zero. Finally, we need to check that Jn{^'i {p}) 
is the reciprocal of a polynomial by calculating its |e| — > oo behaviour. 

Let us illustrate this method to prove that for a + 6 := 2r, c e Z>o, the Fourier integral 
in (1.3) as a function of a — := 2e is equal to the product of gamma functions in (1.3) 
up to a multiplicative function of c and r. For this problem we have 

N 

f{h,...,tr,;{p}) = x{{i-tir{i-i/ur n (i - w^.o^a - w^^r (2.7) 

1=1 l<j<k<N 

and using the functional equation for the gamma function 

TT 



1(^)1(1-^) = -, 



smnz 



we see that the product of gamma functions in (1.3) can be written in the form of (2.6) 
with 

^Me;w)^n' p(^':;:;11c) 

Consider (2.8). Since r and c are integers we see from the recurrence of the gamma 
function 

r{z + i) = Tiz) 

that A^^e; {p}) is the reciprocal of a polynomial, with poles of order j {j = 1, . . . , N — 1) 
at 

|e| = r + 1 + (A^ - j - l)c, r + 2+{N -j - l)c, . . . ,r + c+ {N - j - l)c (2.9) 
and poles of order N at 

|e|=0,...,r (2.10) 
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The order of the polynomial is equal to 

N 
n=l 

where Z„ is the number of poles of order n. Since Z„ = 2c (n = 1, . . . , N — 1) and 1^ — 2r+l 
the order is thus 

cN{N -1) + N{2r + 1) (2.11) 

Now consider the r.h.s. of (2.1) with / given by (2.7). Then (2.1) can be written in the 
form (2.5), where in general k{{p}) may be a polynomial in e. To prove that Aisf{e; {p}) 
is given by (2.8) we need to show three features: 

(a) |?^j({j9})| in (2.5) is no greater than r + c{N — 1). 

(b) For |nj({p})| given by the r.h.s. of (2.9) and (2.10), qj{{p}) is less than or equal to 
the order of the corresponding pole in (2.9) and (2.10). 

(c) For large-|e|, Jn{^', {p}) has an inverse power law decay with exponent (2.11). 

To check (a), we see from (2.7) that the largest power of say ti is r in the product 
over / and c{N — 1) in the product over j and k. It is thus r + c{N — 1) in total, which 
proves (a) for nj{{p}) positive. Since (2.7) is also unchanged by the replacement tj i— > l/tj 
{j — 1,. . . ,N), the same is true for nj{{p}) negative. 

Consider statement (b). With q^^\n) denoting the maximum number of times the 
exponent n occurs in any term of the Laurent expansion of (2.7), statement (b) says 

q^^\n) < N for \n\ < r (2.12a) 

q'-^\n) <N-j for \n\ = r + (j - l)c + z/ (2.126) 

where j = 1, . . . , N — 1 and u = 1, . . . , c. Since there are only N variables g(^)(n) < N 
for any n so (2.12a) is true. To prove (2.12b) consider the explicit formula (2.7) for /. 
Since the first two products can create all powers |n| < r in each variable independently, 
we see (2.12b) is equivalent to saying 

q^^\n)<N-j for \n\ = {j - l)c + u (2.13) 

where q^^\n) denotes the maximum number of times the exponent n occurs in the Laurent 
expansion of 

n (1 - h/mi - tj/t^r 

l<j<k<N 

^(_l)ciV(iV-l)/2-Q^-c(iV+l)+2c,- -Q (l_i,/i.)2c (2.14) 
j=l l<j<k<N 

Writing 

^ = i^...^, j<k (2.15) 

tk tj+2 tk 

we see by expanding the last product in (2.14) that the terms in the Laurent expansion 
of (2.14) are of the form [9] 

N 

^7'=(-^+^+2cj+"j+i-"j (2 16) 
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where rii = un+i = and nj > for each j = 2, . . . , N. Equivalently, setting 

Up+i = c{N - p)p + mp+i, mp+i > -c{N - p)p (2-17) 

for each p— 1, . . . , A?" — 1 we have that all terms in the Laurent expansion of (2.16) are of 
the form 

r/'^'QI''^''-' . . . t^p+^-'^p . . . (2.18) 

We want to determine the maximum number of exponents in (2.18) which can take 
the value (j — l)c+ u. Since (2.14) is a symmetrical function of all the variables we can 
suppose that the k variables t^, tN-i, ■ ■ ■ , tN+i-k have exponent (j — l)c + u. Then from 
(2.18) we require 

mp+i = -((j - l)c + i/)(7V - p), p^N-l,...,N-k (2.19) 
Combining (2.19) with the inequality in (2.16) gives 

(j -l)c + p <cp< c{N - k) (2.20) 
which, since 1 < i' < c implies 

k<N-j (2.21) 

The inequality (2.21) is precisely the statement (2.13) with n = {j — l)c + u. Also, since 
(2.14) is unchanged by replacing all variables by their reciprocals, (2.21) also establishes 
(2.13) for the remaining case n = —{{j — l)c + p). 
To check (c) we note from (2.7) and (1.3) that 



Me; {p}) = (-i)'-M(iv-i)/2 1^1 dti t-i+-'-(^-i)^(i - ti) 

\l=l 



2r 



X n itk-tjf' (2.22) 

l<j<k<N 

The change of variables 

ti = e~^^ then si i— > si/e 
gives the large-|e| asymptotic behaviour 

Me;{p}) - (2.23) 



hN{{p}) := i-ir^<^^M/2 (f[ r ds, sTeA n - s,f 

\/=l"^° / l<j<k<N 



with 

/ N \ 

l<j<k<N 

which is precisely the inverse power law decay with exponent (2.11) required by (c). 

By checking (a)-(c) we have shown, by Liouville's theorem that the Fourier integral in 
(1.3) as a function of a — 6 := 2e is evaluated by the product of gamma functions in (1.3), 
up to multiplicative terms independent of e. 

Remark: The multiplicative function of N,c and r = {a + b)/2 in (1.3) undetermined 
by the above can readily be calculated (see the final paragraphs of section 3.1 for the 
method). We have thus provided a new proof of Morris's integral and consequently of the 
Selberg integral. 
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2.4 Notation 

In the remainder of this paper we will consider the generalization of (1.3) 
Dp{Ni,...,Np;No;a,b, A) 



2mx[°'^ la+b 



Va=l j=l"'-V2 

7 1 — X / ^ / 3 7 • • • J OL 



(2.24) 

where ipo is given by (1.2). We notice that 

Di{N; 0; a, 6, c - 1) = Do{N; a, b, c) = D{N; a, b, c) (2.25) 
where D{N; a, b, c) is the integral in (1.3). 
3. THE CASE p=l 

In the case p — l,a — b — 0, the exact evaluation of (2.24) is known analytically [10]: 



Di(iVi;iVo;0,0,A) 

T{\No + (A + l)A^i + l)r(iVi + 1) 



r(i + A)^or(2 + A)^i (i + S 



(3.1) 



where {a)n :— a{a + 1) . . . (a + n — 1). 

For general a and b the method used in [10] to prove (3.1) docs not appear to be 
applicable. We thus resorted to the numerical approach outlined in Section 2.2. 

3.1 The case A = 1 

With A = IjA^o = 1 and 2, and various values of a and b, we found by sequentially 
increasing Ni that our data fitted the following forms 

Afi-l 

Di(iVi;l; 1,1,1) = 2 n (^■ + l)(2j + 3) 

j=Q 

Ni-l 

Di(iVi;l;2,l,l) = 3 [] (j + l)(2j+4) 

j=0 

D,(iV,;l;2,2,l) ^ ^ ^ 0" + Dl^." + 4)(2, + 5) 
A(7V,;2;1,1,1) = 6 n (j + l)(2j + 4) 

Ni-l 

L>i(7V,;2;2,l,l) = 12 JI {j + l){2j + 5) 

3=0 



j=0 



2j + 4 
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valid for iVi > (when A^i = the products are taken as unity). 

To help fit this data into an analytic form we note from Morris's integral (1.3) in the 
case c — 2 and a e Z>q and h e Z+ that 

^ ^ /i (2j + 2)(2j + 3) . . . (2j + 6) ^ ' 

We see that the above data fits a similar form: 

Di(iVi;iVo;a,6,l) 

f(N h\ "tt' + + « + iVo + l)(2j + g + iVo + 2) . . . (2j + g + 6 + iVp) 
/I o,a,oJll (2^- + Aro + 2)(2j + iVo + 3)...(2j + iVo + 6) 

^fr^ (j + l)r(2j + a + h + No + l)r(2(j + 1) + iVo) 
- ^^"^"'"'^^ M r(2, + g + iVo + l)r(2, + 6 + iVo + l) ('-'^ 

To evaluate f{No, a, b) we set A^i = in (3.3) (the product over j is then taken as unity) 
to obtain 

/(TVo, a, h) = Di(0; N^- a, b, 2) = D{No; a, b, 1) (3.4) 

where D{Nq; a, b, 1) is given by (1.3). 

The resulting conjecture for Di can be proved using the method of Section 2.3. Thus 
we consider Di{Ni; Nq; g, b, 1) as a function of a — 6 := 2e. We have the following result. 

Theorem 1 

Suppose r e Z>o and let 

f{ti, . . . , t]Vo; si, . . . , SNi;r) 

No Ni 

11(1 -wi-i^r 11(1 -«o^(i-iMr 

X n {i-tk/tj){i-tj/h) n (i-ws.)'(i- vsfe)' 

A'o A'l 

xnn(i-w^.)(i-w«fe) (3-5) 

j=i k=i 

Then for Re(e) large enough so that the l.h.s. converges, 

n dtitr'+' \ m dsisf'+'\f{h,...,tMo]s,,...,SM,]r) 

r(e-r-l) ^T^' r(e-r-2l-No) 

where A(A^o, -^i, '^) is independent of e. 
Proof 

Since r e Z>o we see that the r.h.s. of (3.6) is the reciprocal of a polynomial in e. 
Furthermore this polynomial is naturally factored as two polynomials, one for each of 
the products. The first reciprocal polynomial factor is precisely (2.8) with c = 1 and 
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N = Nq. This factor therefore has poles at (2.9) and (2.10) with the order given therein. 
The second factor has poles of order j {j = 1, . . . , Ni — 1) at integers e satisfying 

r + No + 2{Ni - j) < |e| < r + + 2{N^ + l (3.7) 

and poles of order A^i at 

|e| = 0,l,...,r + Aro (3.8) 

The poles at (3.7) and (3.8) occur at the same points as the poles of the first factor, 
plus some additional points. Thus if we let Q'^^'^'^'^\n) denote the order of the pole of 
(3.5) at n, {n & Z), we see that 

Q(^°'^i)(n) = A^o + Ni, \n\ < r (3.9a) 

Q(iVo,ivO(^) = AT, + - j,, |n| = r + ji (3.96) 

where ji = 1, . . . , A^o, 

Q(^o,iVi)(^) = AT, _ \n\=r + No + 2{j2 - 1) + 1/2 (3.9c) 
where j2 — 1, . . . , Ni — 1, 1/2 — 1, 2, and 

g(^°'^i)(n) = 0, otherwise. (3.9d) 
From (3.9) we see that the reciprocal of the r.h.s. of (3.6) is a polynomial of order 

7Vo-l A^i-l 

{2r + l){No + N^) + 2 ^ (j + N^) + A ^ j 

j=0 j=l 

= (2r + l){No + Ni) + 2NoNi + {Nq - l)No + 2{N, - 1)N, (3.10) 

Let g(^"'^i^(n) denote the maximum number of times the exponent n can occur in a 
term of the Laurent expansion of (3.5). In accordance with the method of Section 2.3 we 
want to show that 

^(iVo,iVi)(^) < g(^'''^i)(n) (3.11) 
First, for \n\ < r there is nothing to prove as (3.11) reads 

g(^°'^i)(n) <No + Ni 

which is true by definition. Next consider the statement (3.11) for \n\ and Q^^°'^^\n) given 
by (3.9b). Since the products over / in (3.5) give all exponents n of ti, . . . , tAr„, Si, . . . , SjVi 
with \n\ < r in each variable independently, and (3.5) is unchanged by replacing each 
variable by its reciprocal, we see that in these cases (3.11) is equivalent to proving 

g(^0'^i)(n) < A^i + A^o - Ji, n = j, (3.12a) 

and 

giNo,N,)^n) <Ni- j2, n^No + 2{j2 - 1) + U2 (3.126) 

where q^^°'^'^\n) denotes the maximum number of times the exponent n can occur in the 
Laurent expansion of 

n {^-tk/tj){i-tj/tk) n i^-sk/sjf{i-sj/skf 

l<j<k<No l<j<k<Ni 
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No Ni 

xIlEi^-Sk/tj){l-tj/s,) (3.13) 

j=i k=i 

Using a confluent form of of the Vandermonde determinant expansion, we have previ- 
ously shown [11] that (3.13) is equal to (up to an unimportant ± sign) 

m Ni 

E <n E <Q) n ^>^(^-^)-^-^-^(p(2j) - p(2, - 1)) 

P Q=l j=l 

P{2l)>P(2l-l) 

No 

^ -Q ^P(2iVi+fe)+Q(fc)-JVo-iVi-l 
fe=l 

where P is a permutation of {1,2,..., 2Ni + Nq} with parity e{P) and Q is a permuation 
of {1, ... , Nq} with parity e(Q). 

Let ji, k E {1, . . . , Nq} and j e {1, . . . , A^i}. We see from (3.14) that an exponent of 
ji in the variable requires 

P{2Ni + k) ^ji + No + Ni + l- Q{k) (3.15a) 

and thus 

P(2iVi + A;) e {ji + iVi + 1, ji + iVi + 2, . . . , min(2iVi + Nq.j^ + A^i + Nq)} (3.156) 

For an exponent of j'l in the variable Sj we see from (3.14) that we require 

P(2j) + P(2j - 1) = 2N, + TVo + 1 + Ji, i"(2j) > P(2j - 1). (3.16) 

Consider the case ji < A^i. Then 

min(2Ari + Nq, ji + Ni + Nq) = ji + Ni + Nq 

and so the maximum number of solutions of (3.15a) is A^o- When (3.15a) has this maxi- 
mum number of solutions, (3.16) has solutions for 

(P(2j),P(2j-l)) = (ii + l + iVo + iVi,iVi), (ji + 2 + iVo + iVi,iVi-l),...,(2iVi + iVo,ji + l), 

thus giving a maximum of A^i — ji exponents ji to the variable tj, and thus a total of 
Nq + Ni — ji exponents ji in (3.14). Since decreasing the number of solutions of (3.15a) 
by say Oi can give no more than oi new solutions to (3.16), we thus have shown that for 
ji < Ni 

< ^1 + ^0 - ji (3.17) 
In the cases that ji > Ni (which requires Nq > Ni) 

min(27Vi + Nq, ji + Ni + Nq) = 2Ni + Nq 

and so the maximum number of solutions of (3.15a) is Ni + Nq — ji. When (3.15a) has 
this maximum number of solutions, there are no solutions to (3.16). Arguing as in the 
sentence including (3.17), we conclude that (3.17) also holds for ji > Ni and thus (3.12a) 
is true. 
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For an exponent of n = iVo + 2(j2 — 1) + in the variable sj, (3.14) gives that we 
require (3.16) with ji replaced by n. The maximum number of solutions occurs with 

(P(2j), P{2j-l)) = (27Vi+7Vo, 1+n), (2iVi+7Vo-l, 2+n), . . . , {N,+No+j2+l, N^-j2+n) 

(note that Ni — j2 -\- n = Ni + Nq -\- j2 + 1/2 — 2) and is thus equal to A^i — j2. For an 
exponent of n = No + 2(j2 — 1) + z/2 in the variable tk, (3.14) gives that we require (3.15a) 
with ji replaced by n. When (3.16) has its maximum number of solutions we see that 
(3.15a) doesn't have any solutions. Furthermore, by decreasing the number of solutions 
of (3.16) by say ai, we see that the number of solutions of (3.15a) can increase by no more 
than ai (since for (3.15a) to have a solution we require P{2Ni + k) > n + Ni + 1) and so 
we conclude that (3.12b) is true. 

The validity of (3.12) implies that as a function of e, the l.h.s. of (3.6) divided by 
the r.h.s. is bounded in the finite plane. Furthermore, using the method given in the 
paragraph including (2.22) above, it is straightforward to show that for large-|e| the r.h.s. 
of (3.6) has a reciprocal power law decay with exponent (3.10), which is the same as the 
large-|e| behaviour of the l.h.s. (recall the sentence including (3.10)). Hence, by Liouville's 
theorem, both sides of (3.6) are the same functions of e, up to a multiplicative function 
independent of e. This is the required result. 

Let us now specify the dependence on r of the function A{Nq, Ni,r) in (3.6). For 
this purpose we observe that when e = r, the l.h.s. of (3.6) is independent of r. More 
explicitly. 

No Ni 

CTllf^l[sU{-ti,...,-tNo;-si,...,-SM,;r) 

1=1 1=1 

No Ni 

= CT n(l + 11(1 + si^fi-h, -s^, . . . , -s^,;0) 

1=1 1=1 
— CT/(— ii, . . . , —tNo', —si, . . . , — SiVi; 0) 

where the last line follows from the second last line after noting / is a homogeneous 
function of order 0. For the r.h.s. of (3.6) to have this property, we require 

No-l Ni-1 

A{No, A^i, r) = B{No, N^) [] r(2r + 1 + /) [] r(2r + 1 + 21 + Nq) (3.18) 

1=0 1=0 

The remaining unknown function B{No,Ni) can be specified immediately from the 
anlaytic result (3.1), or alternatively by using the general relationship 

Di{Ni; No; A, A, A) = Di{Ni, No + 1; 0, 0, A) (3.19) 

together with (2.25). We find 

iVo-l ATi-l 

B{No, N,) = n r(Z + 2) n {j + l)r(2(i + 1) + No) (3.20) 

Substituting (3.20) in (3.18), then substituting the resulting expression in the r.h.s. of 
(3.6), we obtain the exact evaluation of the integral in (3.6). This exact evaluation agrees 
with the conjecture (3.3), and thus proves the conjecture. 
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3.2 General A 

Using Morris's integral (1.3), the result (3.4) for A = 1, and the analytic result (3.1) as 
guides, we conjecture that for general A 

D,{NuNo;a,b, A) 

D(N h \) rr' + ^^^^^^ + i)j + g + 6 + xnq + i)r((A + i)u + 1) + xnq) 

^ °' ' ' ^ r(l + A)r((A + l)j + a + AAro + l)r((A + l)j + 6 + AiVo + l) 

(3.21) 

where D{Nq; a,b, X) is given by (1-3). As well as being consistent with theorems used 
in its formulation, (3.21) satisfies the general relationship (3.19), and the consistency 
of the large-|e| behaviour of both sides (recall Section 2.3) can be checked as can the 
independence on a when 6 = (recall the paragraph above (3.18)). 

3.3 Normalization of the harmonic well wavefunction for p = 1 

Prom the conjecture (3.21), it is possible to deduce the value of the integral 

G,{N,,...,Np-No;X) 

(fi ri r dxf] f n r dy,] i#^({xf }.^-...^. ,{%-'},'=i,...,ivo)r 

(3.22) 

where i/jq^^ is given by (1.4), in the case p — 1. 

By setting a = 6 in (2.24) and changing variables xf^ i-^ x^^^ /27ra, yi i— > yi/27ra, we 
see that for a — > oo 

Di(Ni;No;a,a,X) 

~(^) (^) 2^«(^i+^"Gi(iVi;iVo;A) (3.23) 

On the other hand, with a = 6 it is straightforward to obtain the large-a behaviour of the 
conjectured evaluation (3.21) of Di by using Stirling's formula. Comparison with (3.23) 
then gives 

G.liV, iV.; A, ^ (2.,--../^ "fl' ^<V^M 'jj' mX^l)iMlHXN.) ^^^^ 

as the conjectured evaluation of Gi (in the case A = 1, since we have proved (3.21), we 
also have proved (3.24)). 

4. THE GENERAL p cases 

4.1 The case p = 2 

With p = 2 we have the general relations 

L>2(0, N2; 0; a, b, A) = D{N2; a,b,X+ 1) (4.1) 
where the r.h.s. is given by Morris's integral (2.3), and 

^2(1, N2; 0; 0, 0, A) = D,{N2; 1; 0, OA) (4.2) 
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where the r.h.s. is given by the conjecture (3.21). 

These formulas provide analytic data for the cases A^o = and 1. For Ni — 2 and 3, 
and with A = 1, a = 6 = 0, numerical data was obtained. By sequentially increasing N2, 
the data was seen to fit the forms 

16 ^2 

D2(2,Ar2;0; 0,0,1) = -l[j{2j + l), TV^ > 1 (4.3) 

N2 

D2(3,iV2;0; 0,0,1) = 70nj(2j + 2), TV^ > 2 (4.3) 

j=i 

The results (4.1)-(4.4) suggest that 

N2 

D2(iVi,iV2;0; 0,0,1) = g{N,) l[j{2j - 1 + N,) 

r(2j + i + 7Vi) ' ' ^ ^ 

For general a, 6 and A inspection of (4.5) and use of (4.1) and (1.3) suggest the same 
ansatz used in (3.3) for Di{Ni; Nq; a, b, 1): 

D2{N^, N2; 0; a, b, A) = f{N^, a, b)A{N^, N^, a, b, A) (4.6a) 

where 

A(N N h W V (j + l)r((A + l)j + a + b + N, + l)r((A + l)(j + 1) + iVQ 
A(iV„Ar,a,M)= n r((A + l),+a + iV, + l)r((A + l), + 6 + iV, + l) 

(4.66) 

valid for 

A^2 > A^i - 1 (4.6c) 

The restriction (4.6c) is the key distinguishing feature between (3.3) and (4.6a). The 
function f{Ni,a,b) can be determined from the symmetry relation 

D2{N, AT - 1; 0; a, b, A) ^ D2{N - 1, N; 0; a, b, A) 

which gives the difference equation 

f{k; a, b, X)A{k, k-l;a, b, A) = f{k - 1; a, b, X)A{k - 1, k; a, b, A) 

This difference equation has solution 

j^^AT u ^^ ¥\ A{k-l,k;a.b. X) , 

^^^-^"'^'^^^n ^(M-i;M,A) ^'-'^ 

where we have used the fact that 

/(0;a,6,A) = l 

which follows by choosing A^i = iV2 = in (4.6a). Substituting (4.7) in (4.6a) gives the 
conjectured evaluation of D2{Ni, N2; 0; a, b, A). 
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4.2 The general case 

Guided by (4.6a) and (3.3), for general p and A^o we conjecture that 

Dp{N^, ...,Np; No; a, b, A) = /p_i(iVi, . . . , iVp_i; Nq; a, 6, X)Ap{N^, . . . , AT^; Nq; a, 6, A) 

(4.8a) 

where 

Ap{Ni,...,Np;No;a,b,X) 

^i^^ (j + l)r((A + l)j + a + 6 + A E?Io + l)r((A + l)(j + 1) + A E^^o N,) 



- n 



^0 r(l + A)r((A +l)j + a + X E?=o Nj + l)r((A + l)j + 6 + A E?=o A^.- + 1) 

(4.8b) 

and Np > Nj — 1 {j — 1, . . . ,p — 1). To calculate /p_i we make the ordering 

Nj>Nj_^ {j = 2,...,p) (4.8c) 

and use the symmetry relation 

Dp{Ni, Nk-2, N-1,N, Nk+i, ...,Np; Nq; a, b, A) 

= Dp{Ni, Nk-2, N,N-1, Nk+i, ...,Np; No] a, b, A) (4.9) 
for A; = 2, . . . ,p. Prom (4.9) and the initial condition 

Dp{0, . . . , 0; iVo; a, b, A) = D{No] a, b, A) 

where D{No; a, b, A) is given by (2.3), we obtain the recurrence equations 

fk-i{Ni,...,Nk-i;No;a,b,\) 

= Ak-i{Ni, Nk_i; No; a, b, A)/fc_2(iVi, . . . , N^.^; N^; a, b, A) (4.10a) 

where 

A-i(A^, . . . , A^-i; ^o; a, b, A) := U ' ' ' ' ''''^ ^ (4-106) 

fJl Ak{Ni, Nk-2, 3,3 - I; No; a, b, A) 

and 

fo{No;a,b,\) = D{No;a,b,\) (4.10c) 

Taken in the order k = p,p — 1, . . . ,2 these equations explicitly determine fp-i and thus 
Dp. 

For example, with p = 3 we obtain 

D,{N^,N2,N,;No;a,b, A) 

= D{No; a, b, X)As{Ni, N2, Ns; No; a, b, A [[ . . — — - 

j=i A3{Ni,j,j - l;No;a,b,\) 

^ A,{k -l,k-l,k;No; a, b, A) t^-^ A^jk - 1, j - 1, j; Np; a, b, X)A,{k,j,j - 1; Np; a, b; A) 
feii Mk -l,k,k-l; No; a, b, A) /i, As{k - 1,3,3 - 1; ^0; a, b, X)As{k,3 - 1,3; a, b, A) 

(4.11a) 
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where 

Np>Np_,-l, p=2,3 (4.116) 

Note that this agrees with (4.6a) when Ni = Nq = 0. Also, we have made the exact 
numerical evaluations 

L'3(l,2,2;0;0,0,l) = 720 and £13(2, 2, 2; 0; 0, 0, 1) = 10, 080 

and found agreement with (4.11). 

4.3 Normalization of the harmonic well wavefunction in the general case 

Analogous to (4.6a), for the integral (3.22) we conjecture 

Gp{Ni, ...,Np; No; A) = gp-i{N,, . . . , Np_,; No; X)Bp{N,, . . . , AT^; A^o; A) (4.12a) 
where 

B,{N^, ...,N,; No; A) (27r)^^/^ J] fWTTT^^^^ + + 1) + A E ^.0 (4-12&) 

With the ordering (4.8c), and assuming the analogue of the symmetry relation (4.9), the 
conjecture (4.11a) gives the recurrence equations 

gk-i(Ni, Nk-i; No; A) = Bk-i(Ni, Nk-i; Nq; X)gk-2(Ni, ■ ■ ■ , A^fe-2; No; A) (4.13a) 
where 

R (M M ■ V ■ W- • • • ' ^k-2j - I J; No; A) (A^'ih\ 

Bk-i[Ni,...,Nk-i,No,A) .- [[ -B-jT-r 77 :r- .^ (4.136) 

j=i Bk[Ni, Nk-2, J,J -i-; No; A) 

and 

goiNo; A) = (2.)-o/^ ^'^r(l+V'^^ ^'-'"'^ 
which when taken in the order A; = p,p — 1, . . . , 2 explicitly determine gp-i and thus Gp. 
5. SUMMARY 

The objective of this paper has been to provide the exact evaluation of the trigonometric 
integral (2.24). This integral is a generalization of Morris's integral (1.3) (which is equiv- 
alent to Selberg's integral (1-1)), and includes as a special case the normahzation of the 
multicomponent wavefunction (1.2). We have been partially successful in this task in that 
(4.8) and (4.10) provide the conjectured exact evaluation of (2.24) expressed in the forms 
of recurrence equations. Moreover (4.8) provides a conjecture for a specific functional 
property of the integrals (2.24), from which their exact evaluation follows. 

In Section 2.3 we have also provided a new proof of Morris's integral, which was used 
in Section 3.1 to prove the conjectured evaluation of Di{Ni; Nq; a, 6, 1). However we have 
not been successful in providing a proof in the general case, which we leave as an open 
problem. 
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